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The radiation energy component of the Hubble function 
and a ACDM cosmological simulation 
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ABSTRACT 

We study some effects the inclusion of the radiation energy component in the universe, 
S7 r , can have on several quantities of interest for the large-scale structure of the universe 
in a ACDM cosmological simulation; started at a very high redshift (z = 500). In 
particular we compute the power spectrum density, the halo mass function, and the 
concentration-mass relation for haloes. We find that tt Y has an important contribution 
in the long-term nonlinear evolution of structures in the universe. For instance, a 
lower matter density power, by ^50%, in all scales is obtained when compared with 
a simulation without the radiation term. Also, haloes formed with the tt r taken into 
account are ~20% less concentrated than when not included in the Hubble function. 
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1 INTRODUCTION 

In the current cosmological paradigm, structures form in the 
Universe by amplification of primordial fluctuations driven 
by a gravitational instability in the expanding universe 
(e.g. Peacock 1999, Weinberg 2008, Mo, van den Bosh and 
White 2010). The growth of the instability can be studied 
analytically and to some extent an exploration of the weak 
non-linear regime may be done. However, the full non-linear 
evolution is essentially studied by means of numerical sim- 
ulations. 

The formation of non-linear structures in the universe 
has been studied by means of cosmological simulations, 
since some of the first works of, for example, Mellot et al. 
(1983) and Davis et al. (1985). Recent simulations have 
reached a high degree of complexity, both those includ- 
ing only dark matter (e.g. Springel et al. 2005, Diemand 
et al. 2008, Boylan-Kolchin et al. 2009, Heitmann et al. 2010, 
Klypin et al. 2011, Prada et al. 2012) and those with bary- 
onic physics included (e.g. Wise et al. 2012, Bird et al. 2013). 

Cosmological simulations considering physics outside 
the standard ACDM cosmology, such as including warm 
dark matter (e.g. Colm et al. 2000) and different kinds of 
equations for the dark matter and energy (e.g. Klypin et 
al. 2003, Linder & Jenkins 2003, Dolag et al. 2004, Grossi & 
Springel 2009, Rocha et al. 2012) have also been increasing 
in complexity and in the physics explored to understand the 
universe. Cosmological simulations with modified dynamics 
have also been done (e.g. Angus & Diaferio 2011). 

In all of the simulations up to now, to our knowledge, 
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there has been no consideration of the radiation energy den- 
sity contribution to the equation of motion of particles in 
a cosmological setting. This is in part understandable since 
the current cosmic background temperature of photons is 
T 7 o = 2.725 K (e.g. Weinberg 2008), leading to an energy 
density of p 7 o = aeT^ = 4.64 x 10 - 34 gcm -3 , and a den- 
sity parameter of photons Q 7 = p 7 o//>oc = 2.47 x 10 _5 /i -2 ; a 
value much smaller than the current matter density param- 
eter Q m « 0.3 and the vacuum energy density parameter 
Qa ~ 0.7. Even if we consider that an additional contri- 
bution to the radiation component of the universe comes 
from the neutrinos from the era of e~ + e + pair annihila- 
tion the situation does not change by much; the total total 
density of radiation (assuming massless neutrinos) becomes 
p r0 = [l + 3 (7/8) (4/ll) 4/3 j p l0 = 7.80 x 10~ 34 gem" 3 lead- 
ing to a radiation parameter Q r = p r o/po c = 4.15 x 10~ 5 h~ 2 . 
Nonetheless, the effect of the radiation energy density be- 
comes more important towards higher redshift s. 

There are several problems, such as the mass function 
at high redshifts (e.g. Reed et al. 2007, Lukic et al. 2007), 
that demand a treatment as accurate as possible of the evo- 
lution of structures in the universe. For example, Reed et 
al. start some of their simulations at z ~ 300, and Lukic 
et al. run simulations going to as high as z — 500 in their 
study. At those starting redshifts for the simulations the ra- 
diation energy density is not negligible, due to its (1 + z) A 
dependence. Also, due to the nonlinear way matter clusters, 
the effect of changing at high redshift the rate of expansion 
of the universe by including a radiation energy term can be 
significant in structures we see today. 

In this Letter we present results of two cosmological 
simulations done within the standard cosmological scenario, 
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but one including the radiation energy density term in the 
equation of motion of dark matter particles. We quantify dif- 
ferences between both cases in regard to the matter power 
spectrum, the mass function and the concentration-mass re- 
lation for halos. Other properties of the clustering of dark 
matter or haloes themselves are not considered here, nor a 
detailed study of each part is considered; such work is post- 
poned for future communications. The objective is to point 
out the need to include the radiation energy term in cos- 
mological simulations, specially those starting at high red- 
shift, for better consistency with the theoretical framework 
of standard cosmology. 

The outline of this work is as follows. In §[2] we describe 
the model used and describe some numerical matters. In 
§ [3] we show some of the results of both of our simulations. 
Finally, in §2] we provide a summary and final comments on 
this work. 



2 MODELS AND NUMERICAL METHODS 

In cosmological simulations the expansion of the universe 
has to be considered (e.g. Hockney & Eastwood 1981). In an 
iV-body simulation with periodic boundary conditions the 
equation of motion of particle i is (e.g. Bertschinger 1998, 
Springel et al. 2001): 
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where the summation goes over all periodic images of the 
particles j, and a is the scale factor of the Universe. The 
evolution of the latter follows from Friedman equation, 
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with H = a I a, Q m the current epoch matter energy density 
parameter, Q r the present day radiation energy density pa- 
rameter, and £1a the vacuum energy density parameter. A 
flat universe has been assumed in the preceding equations. 
Solving the coupled set of equations JT]) and (|2]), along the 
corresponding Poisson's equation, determines the dynamics 
of the iV-body simulation of the Universe. 

We performed two cosmological simulations us- 
ing the publicly available parallel Tree-code Gadget2 
(Springel 2005). This code uses the Hubble function ([2]), 
including a curvature term, but does not include the radia- 
tion component. In GADGET2 such function is required for 
computing the time steps in the advancement of the mo- 
tion of particles, thru drift and kick factors. We essentially 
modified subprograms driftfac.c and timestep.c of the 
Gadget2 code in order for the simulation to account for 
the Q r contribution. We will denote by ACDM, as is cus- 
tomary, the standard cosmological simulation with Q r = 0, 
and with ArCDM the one including the radiation term in 
the Hubble function. 

Our two simulations take as cosmological parameters 
those of the mean values of the Wmap7 results (Komatsu 
et al. 2011), where the matter density Q m = 0.275, spec- 
tral index n s — 0.968, mass fluctuation ag = 0.816 and the 
Hubble parameter h = 0.702, and we take the vacuum pa- 
rameter as = 1 — H m — Q r - The value of Q r used is that 
indicated in Section 



Initial conditions were generated, using a 2nd-order La- 
grangian perturbation code (Crocce, Pueblas & Scoccimarro 
2006), at a redshift of z = 500. The initial linear power spec- 
trum density is calculated using the transfer function from 
the cosmic microwave background code CAMB (Lewis, Challi- 
nor & Lasenby 2000), normalized to the above as value at 
z = 0. The spectrum is evolved back in time to z = 500, 
using the linear growth factor D + J3} given by (e.g. Car- 
roll et al. 1992, Mo et al. 2010), 
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to generate our initial conditions. Computing J3} at z = 500 
for the ACDM case gives 1.93 x 10" 3 while for the ArCDM 
one 1.53 x 10 -3 ; the radiation energy reduces the growth 
factor by w 20% at that redshift. Such difference will be 
reflected in the displacements and peculiar velocities of par- 
ticles at the initial condition. This approach to modifying 
the initial conditions can only be considered approximate, 
but serves our purpose of elucidating differences when in- 
cluding or not the radiation term in (pQ). 

Each simulation box has a comoving length of L — 
100/i -1 Mpc with N p — 512 3 dark matter particles, leading 
to each particle having a mass of ra p = 5.5xl0 8 /i _1 M©. The 
smallest halo we are able to resolve with some confidence has 
a mass of M = 100ra p «6x 10 10 /i _1 M©. The gravitational 
(Plummer equivalent) softening length was kept at the fixed 
value of e = 5/i -1 kpc in comoving coordinates. Halos were 
identified with the AHF public code (Gill, Knebe & Gib- 
son 2004 and Knollmann & Knebe 2009). 



3 RESULTS 

Qualitative structure. The large scale structure in our box 
appears rather similar for both kind of simulations, however 
important differences appear at smaller scales. In Figure [1] 
we show snapshots at different times [z — 2, 1 and z = 0, 
from left to right) of the distribution of dark matter particles 
around the most massive halo (M = 4.5 x 10 14 h~ x M©) in 
our ACDM cosmology (top) and the same region for that of 
the ArCDM (the halo has M = 3.5 x 10 14 h' 1 M ). At z = 
we find for the ACDM cosmology a total of 47,728 halos, 
while for the ArCDM a total of 46,473 is found. 

Power spectrum. In Figure [3] we show the power spec- 
trum density P(k) computed at different redshifts from z — b 
to 2 = 0, for both kinds of cosmological evolution considered 
in this work. It is readily noticeable that more power den- 
sity > 50%, in more than 2 decades in k, is deposited by 
the ACDM model than the ArCDM. The difference tends 
to increase at higher redshifts as shown in the bottom of 
Figure [3] It is worth noticing that the minimum discrep- 
ancy in power density is « 50% even at the smallest scale of 
our simulations. Not including the effect of D+ in the initial 
conditions, just the Q r in the Hubble function, leads to a 
discrepancy of w 15% at z = 0. 

Mass function. The effect of the Q r term in (pQ) can also 
be seen in the halo mass function, F(M) — N^/V 'AlogM 
with A^h the number of halos in a log-bin of mass M and vol- 
ume V] we computed F(M) as in Lukic et al. (2007). The 
mass function of our haloes for different redshifts is shown 
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Figure 1. Projected mass distribution of particles at z = 2, 1 and z = 0, from left to right, in a slice of a ACDM (top) and ArCDM 
(bottom) cosmology. The length of each slice is 10/i -1 Mpc. Points are colored according to local number density. Differences after 
including the Q r term are noticeable in the snapshots; for instance, less substructures appear to be present in the ArCDM case. 




Figure 2. Projected distribution of particles of the three most massive halos in the ACDM cosmology (top) at z = 0, and the corre- 
sponding haloes in the ArCDM one (bottom). The mass of the halos decreases from left to right. Box size is 3/i -1 Mpc on each side. 
Shown are only 10 per cent of the total number of particles, and colored according to their local density as in Figure [T] 
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k [h/Mpc] 

Figure 3. Power spectrum of the simulation with (dashed lines) 
and without (solid lines) the Q T term. At the bottom the frac- 
tional change, in percent, 5P = (Pacdm - -PatCDm)/ ^acdMj is 
shown. Not including f2 r leads to a higher power density at all 
scales; at least to the resolution limit of this study. 
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Figure 4. Mass function, F(M), for halos in a ACDM (solid lines) 
and ArCDM (dashed lines) cosmology. Redshifts corresponding to 
different epochs are indicated in the plot. A lower number of halos 
at higher redshifts are formed in an ArCDM cosmology. 
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Figure 5. Mean concentration-mass relation, c(M), for halos 
found in both of our cosmological runs. Halos in the ArCDM 
(solid lines) cosmology tend to be systematically lower by « 20% 
than in the standard ACDM (dashed lines) at z = 0. 

file (Navarro, Frenk & White 1997) and concentrations are 
computed by the AHF code. We computed the mean c(M) 
relation as in Kwan et al. (2012). In Figure [5] we show the 
mean c-M relation found. Halos at z — formed in the 
ArCDM cosmology tend to be « 20% less concentrated than 
in the standard ACDM, for most of the mass range of our 
halos; at higher masses (~ 10 14 /i _1 M©) the situation is not 
that clear. The ACDM simulations yields results which are 
consistent with those, for example, of Duffy et al. (2008). 
A simple fit to the ACDM results at z = yields c(M) = 
5.2(M vir /Mi 4 )" - 81 , with M i4 = 10 14 /i _1 M©, while for the 
ArCDM halos the coefficient in former expression is 4.7 with 
essentially the same slope. At higher redshifts the slope 
tends to flatten, as found in other works. 

Substructure. The impact of the radiation term also is 
noticed in the number of subhaloes detected by the AHF 
CODE. In the ACDM model the number of subhaloes for 
the six most massive halos are {94,47,64,36,54,51} while 
for the ArCDM we obtained {55, 35, 14, 18, 19, 24}; i.e. the 
inclusion of Q r tends to reduce the number of subhaloes. In a 
future work we explore the substructure differences in more 
detail. 



in Figure [H As noted, the mass function at z = is some- 
what similar with and without the radiation energy density 
term, but the effect of the latter is stronger toward higher 
redshifts; as was also indicated by the behaviour of the P(k). 
Including the Q r term leads to a lowering of the formation 
of halos at higher redshifts. This can have important conse- 
quences for the demographics of haloes that form galaxies, 
and one may speculate that at much higher redshifts in the 
abundance of primordial haloes that would host the first 
stars in the Universe. 

Haloes concentration-mass relation. As a preliminary 
result on the properties of halos formed under the 
two cosmologies explored here, we computed the mean 
concentration-mass relation, c(M), for the halos found in 
our simulations. Haloes are assumed to follow a NFW pro- 



4 SUMMARY AND FINAL COMMENTS 

We have carried out two numerical experiments on the evo- 
lution of Af-body dark matter cosmological simulation, one 
with the usual neglect of the radiation energy density Q r 
term in the Hubble function and one that includes it. 

Different diagnostics, such as the power density spec- 
trum P(k), the halo mass function, and the c-M relation, 
were used to quantify differences that occur when neglect- 
ing the radiation energy density. Other tests, such as the 
subhalo velocity distribution function are not considered in 
this work, in one part due to lacking the mass resolution 
for an adequate treatment and on other due reduced space 
and nature of this Letter. Future works will address some of 
these topics. 
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Including the radiation term makes the structure for- 
mation of the universe to lag in time in comparison to when 
it is not included. This may be understood, at least in the 
linear regime, from recalling the growing mode D+ ©. In- 
cluding Q r makes the expansion rate larger, hence the per- 
turbation growth is reduced. The same kind of effect occurs 
when, for instance, the Qa is included -for instance- in a 
Einstein-de Sitter universe; both terms Qa and Q r enter the 
Hubble function J2} with the same sign. This enhancement 
of the Hubble drag in comparison with the standard ACDM 
treatment has also the effect of reducing in average the con- 
centration of the halos formed; see Figure [5] 

We have shown that including the radiation energy term 
in a high redshift (z = 500) simulation leads to important 
differences in the structures of the universe than when not 
including it. The power spectrum density "deposited" at all 
scales tend to be lower by ~ 50% in the ArCDM cosmology 
than in the ACDM. Also, the effect of the radiation density 
reflects itself in the mean concentration of halos, by lowering 
it ~ 20% at the current epoch. 

All the results presented in this work point toward the 
necessity to include in simulations the Q r term in the Hub- 
ble function. This is particularly important for questions 
regarding the first structures formed in the universe and 
their evolution. Also it may bear importance in problems 
at galactic scale such as the "missing satellites" (Klypin 
et al. 1999, Bullock 2010) or the "too big to fail" (Boylan- 
Kolchin et al. 2011). Implications of the Q r on such problems 
are the subject of future works, as well as comparisons to 
observations. 
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